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We analyze the time reversal of a multiple- 
input-multiple-output (MIMO) system in a space- 
frequency-selective multi-path fading channel de- 
scribed by the stochastic Schrodinger equation 
with a random potential in the strong-fluctuation 
regime. We prove that in a broadband limit the 
conditions for stable super-resolution are the pack- 
ing condition that the spacing among the N trans- 
mitters and M receivers be more than the coher- 
ence length l c and the consecutive symbols in the 
data-streams are separated by more than the in- 
verse of the bandwidth B~ x and the multiplexing 
condition that the number of the degrees of free- 
dom per unit time at the transmitters (~ NB) be 
much larger than the number of the degrees of free- 
dom (~ MC) per unit time in the ensemble of in- 
tended messages. Here C is the number of sym- 
bols per unit time in the data-streams intended for 
each receiver. When the two conditions are met, all 
receivers receive simultaneously streams of statis- 
tically stable, sharply focused signals intended for 
them, free of fading and interference. This indi- 
cates the rough multiplexing gain of NB in channel 
capacity, with the maximal gain per unit angular 
cross section given by BL d i~ d where L is the dis- 
tance from the transmitters to the receivers. We 
show that under the ideal packing condition time 
reversal can result in a high signal-to-interference 
ratio and low probability of intercept, and hence 
is an effective means for achieving the information 
capacity of multi-path channels in the presence of 
multiple users (receivers). 



Introduction 

Time reversal (TR) of waves is the process of recording 
the signal from a remote source and then retransmitting the 
signal in a time-reversed fashion to refocus on the source 



(see 13] and the references therein). The performance 
of TR depends on, among other factors, the reciprocity 
(or time symmetry) of the propagation channel. One of 
the most striking features of time reversal operation in a 
strongly scattering medium is super-resolution, the coun- 
terintuitive effect of scattering-enhancement of time rever- 
sal resolution 0], Q, 0. It highlights the great potential 
of time reversal in technological applications such as com- 
munications where the ability of steering and pinpointing 
signals is essential for realizing the information carrying 
capacity of a multi-path channel as well as achieving low 
probability of intercept Q, [ipf . 

In order to take full advantage of the super-resolution 
effect in a random medium, one has to first achieve sta- 
tistical stability which can be measured by the signal-to- 
interference ratio (SIR) and the signal-to-sidelobe ratio 
(SSR). Statistical stability and resolution are two closely re- 
lated issues that should be analyzed side-by-side; together, 
they are the measure of performance of TR which depends 
on, but is not guaranteed by, the reciprocity (or time sym- 
metry) of the propagation channel. It has been demon- 
strated experimentally that there are at least two routes to 
achieving statistical stability @ , Q . One route is to use a 
time-reversal array (TRA) of sufficiently large aperture; the 
other is to use a broadband signal (even with one-clement 
TRA of essentially zero aperture). There has been many 
advances in analytical understanding of the former situa- 
tion (see : gj and references therein). In many interesting 
applications of time reversal, however, the aperture of TRA 
is typically small compared to the correlation length of the 
medium, therefore the technological potential of time rever- 
sal hinges more heavily on the second route to statistical 
stability. Compared to the case of large aperture the an- 
alytical understanding of the case of broadband signals in 
time reversal has been so far much less complete with the 
exception of a randomly layered medium Q] . 

In this paper we present the time reversal analysis for the 
MIMO broadband channel whose fc-component is described 
by the stochastic Schrodinger equation 
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in the so called paraxial Markov approximation. Here the 
refractive index fluctuation Xz(-) is a J-correlated-in-z sta- 
tionary random field with a power spectral density <l>(p) 
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FIG. 1: MIMO-TRA 

such that E[x z (x)x z '(x')] = 5(z - z') J $(p)e lp (x - x,) dp 
with E standing for the ensemble average; k is the (dimen- 
sionless) relative wavenumber to the center wavenumber 
ko; the Fresnel number 7 = L z j (JcqL^.) is a dimensionless 
number constituting of the center wavenumber ko and the 
reference scales L z and L x in the longitudinal and trans- 
verse dimensions, respectively, see FIG. 1. The notation o 
in eq. Q means the Stratonovich product (v.s. Ito prod- 
uct). For simplicity of presentation we will assume isotropy, 
i.e. $(k) = <&(|k|),Vk € R d and smoothness of <&. 

The stochastic parabolic wave equation is a funda- 
mental equation for wave propagation in a randomly in- 
homogeneous continuum such as underwater acoustic and 
electromagnetic waves in atmospheric turbulence 0, 
as an approximation to the wave equation with random co- 
efficients when backscattering and depolarization are weak. 
It also models the cross-phase-modulation in nonlinear op- 
tical fibers in the wavelength-division-multiplexing scheme 
|l8| . It has a certain degree of universality and encapsu- 
lates the spatial aspect of wave physics, a primary focus of 
the present work, which is lacking in analysis based on ei- 
ther randomly layered media or random matrix theory |3(J ] , 
0h [i3- Here eq. (JTJ is treated as a model equation for 
continuous random media and is studied rigorously here to 
extract useful insights that may shed light on other types 
of complex media. 

Our goal is to show that for the channel described by 
the stability and super-resolution can be achieved si- 
multaneously when there is sufficiently high number of de- 
grees of freedom at the transmitters. In particular, we 
analyze the multiple-input-multiple-output (MIMO) time- 
reversal communication satisfying the packing condition 
that, roughly speaking, the spacing of the M receivers and 
N elements of the time-reversal array (TRA) is much more 
than the coherence length of the channel, and that the 
consecutive symbols in the T-datum streams are separated 
by more than -B -1 , the inverse of the non-dimensionalized 
frequency bandwidth B{= frequency bandwidth xL 2 /L z ). 



Our main theorem says that in the strong-fluctuation 
regime and the broadband limit the MIMO-TRA sys- 
tem achieves stable super-resolution in the sense that both 
the SIR and SSR tend to infinity and that the signal re- 
ceived by each receiver is focused to within a circle of the 
coherence length £ c when the additional multiplexing con- 
dition is also met, namely NB MT(3 C where (3 C is the 
coherence bandwidth. 

To further understand the meaning and implication of 
the result, we may assume without loss of generality that 
the timing of the data-streams are within one interval of 
the delay spread since the signals separated by much more 
than one delay spread are roughly uncorrelated. Because 
the delay spread 8* ~ /3" 1 (cf. the section entitled From 
SISO to MIMO), Tj3 c is roughly the number of symbols 
per unit time in each data-stream and hence MT(3 C is the 
number of the degrees of freedom per unit time in the en- 
semble of all data-streams while NB is the total number 
of the degrees of freedom per unit time at the TRA. The 
multiplexing condition says that the number of degrees of 
freedom of the intended messages must be smaller than 
the number of degrees of freedom available in the channel. 
The main technical ingredient of our approach is the exact, 
universal low Fresnel number asymptotic obtained for the 
two-frequency mutual coherence function. The calculation 
is tight indicating that the multiplexing condition is sharp. 

The main assumption is the 4-th order sub-Gaussianity 
property (|12|l . The Gaussian-like behavior for 4-th or- 
der correlations is widely believed to occur in the strong- 
fluctuation regime, defined by a 2 — D2L 3> 1,0* = 
D2L 3 ^> 1. We will point out some independent evi- 
dences for this in our calculation. Here L is the (lon- 
gitudinal) distance between the TRA and the receivers 
and D-2 = d^ 1 J |p| 2 <I>(p)<ip is the angular diffusion co- 
efficient (hence a* = \[ThL is the angular spread). In the 
strong-fluctuation regime ,7 1 a* is the spread in the 
so called spatial frequency, er* = V D2L 3 the spatial spread 
and their product 7~ 1 £>2-L 2 the spatial- spread-bandwidth 
product (SSB) which, as we will show, is exactly 7" j3~ . 
By the duality principle for the strong-fluctuation regime, 
proved in the effective aperture is 2tt times the spa- 
tial spread er* (independent of the numerical aperture of 
TRA and hence super-resolution), and its dual quantity 
7L/(fccr*) « 7/a* (the inverse of spatial-frequency spread) 
is the coherence length l c of the forward propagation (as 
well as the time reversal resolution). Hence the ratio 
cr*/£ c equals the spatial-spread-bandwidth product and is 
roughly the number of uncorrelated sub-channels (paths) 
per transverse dimension in the cross section of diameter 
o~ * at the receiver plane, which will place upper bound on 
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the capacity gain per unit angle of the channel (see more 
on this in the Conclusion). 

In what follows, we first formulate the problem and 
develop the essential tool for analyzing TR, the one- 
and two-frequency mutual coherence functions, and then 
carry out the stability and resolution analysis for the 
single-input-single-output (SISO), multiple-output-single- 
output (MISO), single-input-multiplc-output (SIMO) and 
the multiple-input-multiple-output (MIMO) cases. Both 
MISO- and SIMO-TRA systems have been demonstrated 
to be feasible for ocean acoustic communication [2^], Q, 
and the MIMO-TRA system with N > M has been 
shown to work well for ultrasound 0. We will discuss the 
implications of our results on the channel capacity in the 
Conclusion. We have by and large neglected the effect of 
noise in our analysis, assuming that the TRA operates in 
a high signal-to-noise ratio (SNR) situation as is the case 
for the experiments reported in 0- The robustness of 
TR in the presence of noises has been well documented, see 




MIMO-Time reversal 

We extend the time-reversal communication scheme Q 
to the setting with multiple users. Let the M receivers 
located at {L,Yj),j = 1, ...,M first send a pilot signal 
/ e l ~ g(k)dkS(rj — an) to the TV-element TRA located at 
(0, a.j), i = 1, N which then use the time-reversed version 
of the received signals J e l ~ g{k)GL{v 3 a,; k)dk to mod- 
ulate streams of symbols and send them back to the re- 
ceivers. Here is the Green function of eq. JJ| and g 2 (k) 
is the power density at k. As shown in |3j, when the 
TRA has an infinite time-window (see the Conclusion for 
the case of finite time-window), the signal arriving at the 
receiver plane with delay L + t is given by 

T N M 

SM) = EEE TO ^ T <) / e-^-^gik) 
1=1 i=i j=i * 

xG L (r,ai; k)G* L (rj , m; k)dk (2) 

where mj(r;), / = 1, ...,T < 00 are a stream of T symbols 
intended for the j-th receiver transmitted at times t\ < 
t 2 < ... < tt- We assume for simplicity that |mj(r;)| = 1, 
V?", I. We have chosen the time scale such that the speed of 
propagation is unity (thus wavenumber=frequency) . 

We assume that g is a smooth and rapidly decaying func- 
tion with effective support of size Bj. For simplicity we 
take g 2 (k) = exp (— ) ■ The broadband limit may be 
formulated as the double limit 

7->Q, B^oo limB7 = (3) 



so that in the limit g 2 {k) becomes narrowly focused around 
k = 1. The idea underlying the definition is to view the 
broadband limit as a sequence of narrow-bands with in- 
definitely growing center frequency and bandwidth. This 
is particularly well suited to the framework of parabolic 
approximation described by Q). The apparent narrow- 
banding of © is deceptive: the delay-spread-band-width 
product (DSB) turns out to be B(3~ x and is doubly diver- 
gent as B — > co (the broadband limit) and f3 c — > (the 
strong fluctuation regime). Note that since k is the relative 
wavenumber, the product Bj should always be uniformly 
bounded between zero and unity, independent of 7 > 0. In 
the case d = 1 this has the intuitive implication that the 
number Bj3~ l of degrees of freedom at each TRA-element 
is less than or equal to the number 7 _1 /3 ( r 1 (SSB) of mi- 
correlated propagation paths in the medium. 

Packing condition. We assume that the spacing 
within the N TRA-elements and the M receivers be much 
larger than the coherence length £ c and that the separa- 
tion of the successive symbols be much larger than (2B)~ l . 
Though there is no technical limitation on M, N, T, it suf- 
fices to consider the case where all the N TRA-elements 
and all the M receivers are located within one circle of 
diameter er* (implying M, N <C "f~ d f3~ d )i an d all the T- 
datum streams are within one interval of the delay spread 
~ /3" 1 (implying T <C B(5~ x ) since the signals separated by 
much more than one spatial spread er* or one delay spread 
5* are essentially uncorrelated. 

For simplicity, we have assumed that all the receivers lie 
on the plane parallel to the TRA. When this is not the case, 
then the above spacing of antennas refers to the transverse 
separation parallel to the TRA. 

SIR/SSR. Anticipating a singular limit we employ the 
coupling with smooth, compactly supported test func- 
tions. Denote the mean by E(r,t) = 7~ d J6*((x — 
r)/£ c )'ES(x,t)dx where the coupling with the test func- 
tion 9 can be viewed as the averaging induced by measure- 
ment. Denote the variance by V(r,t) = j~ 2d E[f 9*((x — 
r)/£ c )S(x, t)dx] 2 — E 2 {r,t). We have made the test func- 
tion 9 act on the scale of the coherence length £ c , the small- 
est spatial scale of interest (the speckle size) in the present 
context. Different choices of scale would not affect the con- 
clusion of our analysis. 

The primary object of our analysis is 

p(r,t) = E ^^ ) , j = l,...,M,l=l,...,T (4) 

which is the SIR if r = r^, t = t\ and the SSR if |r — r 3 - 1 ^> 
^ c ,Vj (spatial sidelobes) or \t — ti\ ^> B^ 1 , V7 (temporal 
sidelobes) (as V(r,r) « E 2 (r,r) as we will see below). 
In the special case of r = rj and \t — r;| ^> B~ 1 ,Vl, 
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p^ 1 is a measure of intersymbol interference. To show 
stability and resolution, we shall find the precise condi- 
tions under which p — > oo and ES'(r, i) is asymptotically 
J2?=iYljLi m j{ T i)Sji{r,t) where Sjt(r,t) = in the spa- 
tial or temporal sidelobes and 



N 



(5) 



g(fc)E[G L (r,a 4 ;fc)Gl(r J ,a 4 ;fc)]dfc 



is a sum of 5-like functions around r, and r ; = 0,V7. In 
other words, we employ the TRA as a multiplexer to trans- 
mit the M scrambled data-streams to the receivers and we 
hope to turn the medium into a demultiplexer by employing 
the broadband time reversal technique. 

Mutual coherence functions 

A quantity repeatedly appearing in the subsequent anal- 
ysis is the mutual coherence function T z between the Green 
functions at two different wavenumbers k\ = k— 7/3/2, k% = 
k + 7/3/2 



x + r x — r 



fc,/3) = E[G,(x,a;fc- 7 /3/2) 
xG:(r,a';fc + 7 /3/2)]. 



We shall omit writing k, j3, a, a' when no confusion arises. 
Here we have chosen x, r to be the pair of variables of con- 
cern and left out a, a' as parameters. By the reciprocity of 
the Green function, we can choose one variable from {x, a} 
and the other from {r, a'} as the variables of T z and leave 
the others as parameters. 

One-frequency version. When /3 = 0, T z satisfies 

3 k 2 

^r, + ^v x ■ v y r z + — D{ 1Y )v z - o (6) 

oz k 7 

where the structure function of the medium fluctuation 
D(x) is given by D(x) = J $(k) [l - e ik x ] dk > 0, Vx e 
~R d . Eq. {HI is exactly solvable by the Fourier transform in 
x. For a* « 7" 1 and 0% <C 7 _1 we can use the approxi- 
mation 

k 2 -r 2 D(-/y-pzj/k)dz^ D 2 \ky-pz\ 2 dz 
Jo Jo 

to obtain 

r L (x,y;fc,0) (7) 
w J e lp x f (p,y- ^;fc,0)e-^ D2 l' £y+pz l 2dz dp 

= J e 1 P x f (p,y-^;fc,0) e -^ly+P 2 l 2d2 dp 



where y = ya*fc and p = per* It is clear from Q that 
has a Gaussian-tail in y (the difference coordinates) and, by 



rescaling, an effective support 



and hence 



7/"* 



(recall that y is the coordinate on the scale 7 _1 ). 

Two-frequency version. The two-frequency mutual 
coherence function is not exactly solvable except for some 
special cases. Fortunately the asymptotic for 7 <C 1 has a 
universal form and can be calculated exactly. Without loss 
of generality we assume /3 > in what follows. 

Using the so called two-frequency Wigner distributions 
we have proved in 9] that in the limit 7 — > 0, Y z satisfies 
the equation 



dT z i 

dz fc Vy Vx z 



-Do 



ky + 



/3 2 

T z - tL. £ r(8) 



where Do — J $(k)(ik. The key to understanding eq. JHJ 
is the rescaling: 



— , y = yfca*, z = z/L, (3 = (3/(3 c 



(9) 



with /3 C = D 2 1 L 2 which transforms eq. JHJ into the form 



9z 



«Vy ■ v x r 2 



|y+2 x l T > 



I3 2 DqT z 
2a 2 



(10) 



By another change of variables zi = y + /3x/2, Z2 = y — 
/3x/2 eq. (|10(l is then transformed into that of the quantum 
harmonic oscillator and solved exactly. The solution is 
given by 



rL(x,y;fc,/3) 



(11) 



e / dx'dy'e & 



sm 



i d / 2 (/3 1 /2(i +i) ) 

-1^4 cot (V^(l+i)) l/Sx+y ^JV I 2 

Xe 2\/^ ' V I cos (V5(l+i)) I 

xe 2 v^' r (o-,x', ). 

ka* 

Several remarks are in order: (i) The Green function for 
is of the Gaussian form in x, y, consistent with the (sub- 
)Gaussianity assumption; (ii) In the vanishing fluctuation 
limit Dq,D 2 — *■ the free-space two-frequency mutual co- 
herence function is recovered; (iii) The apparent singular 
nature of the limit j3 — > in is deceptive. Indeed, the 
small (3 limit is regular and yields the result obtained from 
eq. (jHJ with /3 = 0; (iv) In the strong-fluctuation regime, 
Dq is typically much smaller than D 2 L 3 3> 1 so the factor 
exp ( — P 2 D /(2a 2 )) is negligible in the strong-fluctuation 
regime. On the other hand, the rapidly decaying factor 
sin~ d / 2 (/3 1 / 2 (1 + i)) is crucial for the stability argument 
below; (v) r^(x, y; k, (3) is slowly varying in x on the scale 
<t* for {3 ~ P c and more rapidly varying in x for /3 ^> f3 c . 
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Fourth order sub-Gaussianity. The strong- 
fluctuation regime a* S> 1,(7* 3> 1 can result from ei- 
ther long distance propagation and/or large medium fluc- 
tuation. It is widely accepted that, in this regime, the 
statistics of the wave fields (for at least lower moments) 
become Gaussian-like resulting in, for d = 2, an expo- 
nential PDF for the intensity 0, 0, 0, 0, Q, 
|25| . The Gaussian statistics follows heuristically from 
Central-Limit-Theorem as the number of uncorrelated sub- 
channels (paths) per transverse dimension in the cross sec- 
tion of diameter cr* increases linearly with the spatial- 
spread-bandwidth product, as explained in the Introduc- 
tion. This is consistent with the experimental finding of 
the saturation of intensity fluctuation with the scintillation 
index approaching unity [itI ]. 

In what follows we shall make the 4-th order sub- 
Gaussianity hypothesis, namely that the fourth moments 
of the Green function at different frequencies {Gl(A;)} can 
be estimated by those of the Gaussian process of the same 
covariance. More specifically, we assume that 

|E [G L (ki) ® G* L (h) ® G L (k 2 ) ® G* L (k 2 )} (12) 
-E [G L (ki) ® GUh)} ® E [G L (k 2 ) ® G* L (k 2 )] I 
< K\E [G L (k!) ® G?i(fe)] ® E [G^(fci) ® Gi(*a)] | 
+K\E[G L (ki)(8)Gl(k2)}®E[Gl(k 1 )®G L (k2)} \ 

for some constant K independent of 7 — > 0, |fci — 1| = 
O(B-f), I A^2 — 1 1 = O(Bj) and all the variables. For a jointly 
Gaussian process, the constant K = 1. Note that, in view 
of the scaling in the two-frequency mutual coherence the 
first term on the RHS of (|12fl is much smaller than the 
second term due to difference in wavenumber for Gi(fe) = 
G* L (-k). 

The sub-Gaussianity assumption will be used to estimate 
the 4-th order correlations of Green functions appearing in 
the calculation for V by the two-frequency mutual coher- 
ence function in the strong-fluctuation regime. 



by 



From SISO to MIMO 



Our first application 
ence functions is the 



of the 
estimate 



mutual 
for the 



coher- 
dclay 



spread. 



Consider the band-limited impulse response 

ik(L-t) 

m(x, t) —J g(k)e -< Gl(x, 0; k)dk. It follows easily using 
the preceding results that the mean delay is L and the 
asymptotic for the delay spread 5*, when B 3> (3 C , is given 



= J (t- L) 2 E\u{x,t)\ 2 dt/ / E|u(x,i)| 2 dt 




r L (x,0;l,/3)/r L (x,0;l,0)^/3 c " 1 

/3=0 

which is slowly varying in x on the scale ct* . As commented 
before it suffices to consider the case with a finite T such 
that \t\ — tt\ ~ Pc > implying the number of symbols in 
each data-stream T <C B/3^7 1 , the DSB. In what follows, 
due to ft < 1 the temporal component of the signals is 
essentially decoupled from the spatial component and de- 
termined by the power distribution g 2 . 

SISO. This case corresponds to N = 1,M = 1. Let 
ai = 0. In the calculation of E(x,t), the expression 



?,r z )(r)= / n 



^)r L (r 1 + ^,y;fc,0)dy 



arises and involves only the one-frequency mutual coher- 
ence. Using with r (x,y) = 6(x + ^)<5(x - ^) and 
making the necessary rescaling of variables we obtain the 
following asymptotic 



?,r i )(r)«G (r,r 1 )/3f 
ri — r 

Jr. 



Go = / d P 0*(p 



(13) 

-|p| 2 /3. (U) 



To derive l|13|) we have used the defining conditions of the 
strong-fluctuation regime. Note that the transfer func- 
tion in l|14fl is Gaussian in p and that Go(r, ri) has a 
Gaussian-tail in |r— ri\/£ c and Go(ri,ri) is bounded away 
from zero and slowly varying in ri on the scale <r*. That 
is, after proper normalization Go(r, ri) behaves like a 8- 
function centered at ri. By (|13l) - (|14f) we obtain the mean 
field asymptotic E(r,t) w for |r — ri| 3> l c (spatial 
sidelobes) or \t — r;| 3> B~ 1 ,Vl (temporal sidelobes) and 

The calculation for V involves the four-point correlation 
of the Green functions at different frequencies. Under the 
sub-Gaussianity condition (|12f> the calculation reduces to 
that of two- frequency mutual coherence functions. 

Using (EJ with T (x, y) = c5(x+ ~-f )<5(x- =f) we obtain 
the asymptotic for the dominant term in the calculation for 
U(x, r) prior to the k- integration 



r L (n,o ; fc,/3) Jv L [ 



Xi + X 2 Xi — x 2 



k,0) 



xr(V)«( 



x 2 -r,, Xl ,x 2 _ ^ 



-)d-d- 

«c 7 7 



(15) 
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with the constant Cp given by 

C h = (2TT) 2d (l + i) d p d / 2 sm- d (^/j(l + t))e~ E ^ L 

_L|v'l 2 -l^i cot (\fd(l+ij) I ^ il +y'| 2 , 

Due to the rapidly decaying factor sin~ d (j\J~0(l + £)) 
the /3-integration of Cs is convergent as B — > oo. Be- 
cause /3 C -C 1, in the (fci, fc 2 ) -integration the power distri- 
bution g{k\)g{k2) and are decoupled after the change 
of variables: (fci,fc 2 ) = (k — [3j/2,k + f3j/2), so we con- 
clude that F(x,t)<2V2^^7 2 /3c d+lBT / C^&fi. Note that 
the variance increases linearly with the number T of sym- 
bols in each data-stream. 

The asymptotic SIR/SSR for the SISO-TRA is given by 
p = (Bp- l T- 1 ). Note that the SIR/SSR is slowly vary- 
ing in the test point r and the receiver location ri on the 
scale of ex* . 

SIMO. Let us turn to the SIMO case with N = 1 ele- 
ment TRA located at ai = 0. 

The mean field calculation is analogous to the SISO case. 
Namely, E(Tj,Ti) ps y/A~TrCo(rj,rj)f3^B^/ and zero in the 
temporal or spatial sidelobe. 

In view of the the remark following (|12|) the variance of 5 
is dominated by the contribution from the diagonal terms 
in the summation over receivers given by 

[jot I n ; / x l+ x 2 x l ~ x 2 X x X 2 

x / d/3r L (r J ,0;fc,/3)r L ( , ;k,(3)d — d — 

J 2777 

re VZkB^TM J Cpdppf, 

because |r» — r,-| ^> l c regardless whether the test point 
is near or away from any receiver. Therefore we have the 
estimate: p = O (B fa 1 M^T' 1 ) . 

MISO. The case corresponds to M = 1. Each term 
in the summation over the N TRA-elements has the same 
asymptotic as that of the SISO case. Hence E(rj,Ti) ps 
y/AnNCoi^j, r j)0c^l an d zero in the spatial or temporal 
sidelobes. 

For the variance calculation, let us first note that the 
correlations of two Green functions starting with two TRA- 
elements located at , &j satisfy eq. JSJ) in the variables 
(aj,aj), by the reciprocity of the time-invariant channel, 
and hence vanish as |a$ — &j\ 3> t c - The variance of the 
signal at r (whether at 17 or away from it) before perform- 
ing the fc-integration is then dominated by the following 



diagonal terms in the summation over receivers 

N . _ 

^[G^n.ajsfcOG^n.a^fca)] / 0*£±— -) 
3=1 J c 

x9(^-^)E [G L (xi,a i5 fc 1 )G2(x 2) a j ; fc 2 )] d—d— 

4 7 7 

« NC (3 2d . 

The fc-integration induces the additional factor of 
V2^B-/ 2 T. Hence V(r,t)<2V2^Kj 2 f3 2d+1 BT J C^d(3 
since |r — Tj\ <C cr*,Vj. We conclude that p = 
O (A_B/3~ 1 T~ 1 ). 

MIMO. The analysis for the MIMO case combines all 
the previous cases. The mean signal has the same asymp- 
totic as that of the MISO case, i.e., linearly proportional 
to BN. The variance of the signal prior to performing the 
fc-integration is dominated by 

MM _ 

E[Gt(r i , aj ;k 1 )G L (v haj ;k 2 )] / 

»j'=i c 

x 0(^-^ )E [G L (xi,a,-; fei)Q(x 2) a,; fc 2 )] d^d— 
t c 7 7 

ps NMCpf3 2d 

and therefore F < TMN2K-/ 2 Bf3 2d+1 J C^dp. 

We collect the above analysis in the following statement. 

Summary. Let the N -element TRA, M receivers and 
the number of symbols T satisfy the packing condition. As- 
sume the J^-th order sub-Gaussianity condition in the 
strong-fluctuation regime and let 1 < a, < 7 _1 , 1 <C <J* <C 

Then in the broadband limit ^) the asymptotic 
SIR/SSR ~ NM~ 1 T~ 1 B[3~ 1 is valid uniformly for all 
Tj,j = 1,...,M, with the constant of proportionality 
2- 1 {2ir)- 1 / 2 K- 1 (f Cpdfy-^Col 2 where C as given by 
\14\ is not zero for 6^0. 

The asymptotic signal at the receiver plane within the dis- 
tance <C6% from the receivers is Ylf—i Ylf m j ( T i)Sji ( x ; 
where Sji(x,t) given by @). 

Conclusion and discussion 

The strong-fluctuation regime constitutes the so called 
space-frequency-selective multi-path fading channels in 
wireless communications [22j. In such a channel, TR has 
the super-resolution given by l c — ^j^JD^L. We have es- 
tablished firmly the packing and multiplexing conditions 
for stable super-resolution for the MIMO-TRA communica- 
tion system under the 4-th order sub-Gaussianity assump- 
tion. The experimental evidence for our result in the case 
of M = 1 has been demonstrated in m . 
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We have argued that statistical stability is crucial for 
multi-receiver TR communications, especially when the 
multiple receivers do not have channel state information, 
as the multiuser interference is essentially indistinguish- 
able from the intended signal, the only difference being 
their statistical properties. The latter is in the mean field 
while the former is primarily in the fluctuating field. Our 
result implies that the time-reversal communication can be 
realized stably in principle with up to M ~ NB/3~ 1 T~ 1 
receivers simultaneously at the rate T(3 C with low proba- 
bility of intercept due to super-resolution. Concerning the 
channel capacity, our result is analogous to the finding in 
01- [2l| . [22! based on the random matrix model- 

ing and theory that the ergodic capacity with complete 
channel state information at the receiver with M receive 
antennas (but not at the N transmit antennas) scales like 
min (M, N) log 2 SNR (per unit frequency) at high SNR. Af- 
ter taking into account the frequency multiplexing gain Q , 
the multi-frequency channel capacity then scales like 
B min (M, N) log 2 SNR. Note, however, this result does 
not include the interference due to noncooperating mul- 
tiuser receivers as we do here. Also, these works consider 
only narrow-band signals for which statistical stability is 
rarely valid in practice and consequently the ergodic ca- 
pacity is an average, not almost sure, quantity. 

In the present set-up with the B-band-limited channel 



state information at the transmitters but not the receivers, 
the multiplexing gain is, up to a logarithmic factor, roughly 
MT(3 C ~ BN , the number of degrees of freedom per unit 
time at TRA (see for more analysis on TR capacity 
in multi-path Rayleigh fading channels). The packing con- 
dition also points to the maximal capacity per unit an- 
gular cross section Bl~ d L d log 2 (SNR) when N reaches 
the saturation point a d /£ d in the angular spread a*. Here 
L d £~ d = a dr y~ d has the physical meaning of the angular 
density of uncorrelated propagation paths in the medium. 

Let us point out several possible extensions of our results. 
First, the case of even broader bandwidth of < lim£?7 < 
1 can easily be treated by partitioning the full passband 
into many sub-bands with their own B and 7 satisfying @ . 
Since the self-averaging takes place in each sub-band and 
the whole process is linear, stable super-resolution is valid 
in the full passband. Second, in the case of a finite time- 
window, the out-put signals, unlike J2J), involve a coupling 
of nearby wavenumbers [l2j| . If the time window is suffi- 
ciently large, 3> then the coupling takes place only be- 
tween wavenumbers of separation much smaller than f3 c and 
our result carries over without major adjustment. Finally, 
our results may also be extended to time- varying channels, 
prevalent in mobile wireless communications, with a low 
spread factor T~ 1 5* -C 1 where T c is the coherence time 
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